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ABSTRACT 



We calculate the nucleosynthesis inside the hot bubble formed in the jittering- 
jets model for core collapse supernovae (CCSNe) explosions, and find the forma- 
tion of several xlO _4 M of r-process elements. In the jittering-jets model fast 
jets launched from a stochastic accretion disk around the newly formed neutron 
star are shocked at several thousands km, and form hot high-pressure bubbles. 
These bubbles merge to form a large bubble that explode the star. In the cur- 
rent study we assume a spherically symmetric homogenous bubble, and follow its 
evolution for about one second during which nuclear reactions take place. The 
jets last for about one second, their velocity is Vj = 0.5c, and their total energy 
is ~ 10 51 erg. We use jets' neutron enrichment independent on time, and fol- 
low the nuclear reactions to the formation of the seed nuclei up to Z < 50, on 
which more neutrons will be absorbed to form the r-process elements. Based on 
the mass of the seed nuclei we find the r-process element mass in our idealized 
model to be several xlO _4 M , which is slightly larger than the value deduced 
from observations. More realistic calculations that relax the assumptions of a 
homogenous bubble and constant jets composition might lead to agreement with 
observations. 



1. INTRODUCTION 



The mechanism for the explosion of core-collapse (CC) supernovae (SNe) is still un- 
known. Most popular are models based on ex plosion driven by neutrinos. Less popular 
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Chugai et al.l l2005t ISmith et al.l 1201 ll ) that found asymmetry in CCSNe suggest that jets 
might indeed play an important role in at least some CCSNe. 



I n neutrino-driven models where jets pla y no role at all (e.g. lBethdll990l ; lNordhaus et al. 



2010l : iBrandt et al.l 1201 ll ; lHanke et al.ll201ll ). the neutrinos are absorbed near the stalled- 
shock at ~ 300 km, and revive the shock. Namely, they eject the material in that region. In 
some theore tical studied the jets were injected at large distances beyond the stalled-shock 



radius (e.g ., iKhokhlov et al.lll999l : iHoflich et al.ll200ll : Maeda fe Nomotoll2003l : ICouch et al. 



20091 . l201ll ) . IHoflich et al.ll200ll injected slow and fast jets at a radius of 1200 km in a helium 
star for about two seconds. They show the possibility of exploding the star with slow jets. 
Their result s show that the po larization in SN199em is consistent with slow jets. In model 
m2rlhot of ICouch et al.l 1201 ll a jet was i njected near the speed o f sound, leading to the 
formation of a hot bubble. In others, like iMacFadyen et al.l (120011 ) . the jets were injected 
much closer to the neutron star (NS), at 50 km. In the simulations of IMacFadyen et al. 
( 1200 ll ) the jets were injected at a muc h later time i n the explosion, and are less relevant to 
our goal of exploding a star with jets. iKohri et al.l (120051 ) propose that disk-wind energy is 
able to revive a stalled shock and help to produce a successful supernova explosion. 

Our jittering-jet model for explosion (Soker 2010; Papish & Soker 2011, hereafter Paper 
1 ) is based on the following points, that differ in several ingredients from the models cited 
above (for more detail see Paper 1). (1) We don't try to revive the stalled shock. To 
the contrary. Our model requires the material near the stalled-shock to fall inward and 
form an accretion disk around the newly born NS or black hole (BH). (2) We conjecture 
that due to stochastic proce sses and the stationary accretion shock instability (SASI; e.g. 
Blondin fc Mezzacappall2007l ) segments of the post-shock accreted gas (inward to the stalled 
shock wave) possess local angular momentum. When they accreted they form and accretion 
disk with rapidly varying axis direction. (3) We assume that the accretion disk launches 
two opposite jets. Due to the rapid change in the disk's axis, the jets can be intermittent 
and their direction rapidly varying. These are termed jittering jets. (4) We show in Paper 
1 that the jets penetrate the infalling gas up to a distance of fewxlOOO km, i.e., beyond the 
stalled-shock. However, beyond fewxlOOO km the jets cannot penetrate the gas any more 
because of their jittering. The jittering jets don't have the time to drill a hole through the 
ambient gas before their direction changes; they are shocked before penetrating through the 
ambient gas. This condition can be met if the jets' axis rapidly changes its direction. This 
process of depositing jets' energy into the ambient medium to prevent further accretion is 
termed the penetrating jet feedback mechanism. (5) The jets deposit their energy inside the 
star via shock waves, and form two hot bubbles, that eventually merge and accelerate the 
rest of the star and lead to the explosion. In section [2] below we use self similar calculations 
to further explore this process. (6) The jets are launched only in the last phase of accretion 
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onto the NS. For the required energy the jets must be launched from the very inner region 
of the accretion disk. 

Nucleosynthesis can occurs i n the exp a nding jets (or disk winds) and in the postshock 
region. Previous studies include ICameronl ( 120011 ) who discussed the nucleosynthesis inside 
jets launched at a velocity of 0.5c from an accretion disk around a rapidly r otating proto NS. 



He su ggested the possibility of creating r-process elements inside the jets. iNishimura et al. 



(120061 ) simulated the r-process nucleosynthesis during a jet powered explosion. In their 
simulation a rotating star with a magnetic field induces a jetlike outflow during the collapse 
which explode the star. Neutrinos play no role in the simulation. Unlike their model, our 
model does not have a large scale rotation of the star, and the jets penetrate farther away 
creating hot bubbles. 

In our jittering-jets explosion model the jets are launched close to the NS where the 
gas is neutron-rich (e.g., iKohri et al.l 120051 ) . In section [3] we examine the implications of this 
on the nuclear reactions (nucleosynthesis) in the inflated hot bubble. The properties of the 
bubbles are as derived in section |2J Our discussion and summary are in section HI 



2. DYNAMICAL EVOLUTION OF THE INFLATED BUBBLES 

In this section we describe an approximate model for the inflat ed bubbles an d their 



dynamical evolution. The analysis here is similar to that conducted by IVolk fc Kwokl (119851 ) 
to study the evolution of the spherical hot bubble in planetary nebulae. During the active 
phase of the jets we derive a self-similar analytical solution to the gas-dynamical equations. 
At later times the solution is numerical. 

The jittering jets form wide bubbles that occupy most of the volume up to the distance 
they have reached; eventually the bubbles merge (see Paper 1). Our basic assumption is 
therefore that the two inflated bubbles merge to form one large bubble. The low density 
high energy volume inside R s is termed hereafter the 'spherical bubble'. If we equate the 
volume of the assumed spherical bubble with the total volume of the wide inflated bubbles, 
the radius of the spherical bubble R s is only slightly smaller that the distance the inflated 
bubbles have reached. This assumption leads to a spherically symmetric flow that allows a 
self similar solution for constant power jets. We also assume that the gas inside the bubble 
is homogeneous, i.e., the composition, density, and pressure are constant inside the bubble. 
The energy of the jets injected into the bubble with a power of Ej forces the bubble to 
expand. The expanding spherical bubble pushes the dense surrounding gas supersonically 
outward, forming a forward shock ahead of this dense shell. The boundary of the dense 
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shell and the bubble is a contact discontinuity, across which the pressure is constant but not 
the density or the composition. This flow structure is drawn schematically in Fig. [TJ The 
mass M s in the shell is the swept-up ambient gas. As the dense shell is thin, we take the 
radius of the forward shock to be equal to the radius of the contact discontinuity (which is 
the radius of the spherical bubble) R s . The pre-shock (up-stream ) ambient density p rofile 
at radius r > is taken to be a power law, with the scaling from I Wilson et al.l (119861 ) and 



Mikami et all (120081 1 (see Paper 1) 
p s {r) = Ar 13 = 1.3 x 10 



10 



100 km 



-2.7 



g cm" 3 , 30 < r < 10 4 km. 



The spherical flow obeys the following conservation equations (jVolk &: Kwoklll985l ) 

<1}L 4nR*p(R s )R s , 



dt 



d_ 

dt 



j t (m s R s ) = 4ttR 2 s PR s , 
(4ttR 3 s P) = Ej - 4ttR 2 s PR s 



(2) 

(3) 
(4) 



where P is the pressure inside the bubble. Equations ([2j) - (j4]) describe the conservation 
of mass, momentum, and energy respectively. The energy inside the bubble includes the 



Jets termination 
shock 



Forward shock 



Stellar core 




Stellar core 



Contact 
discontinuity 

Fig. 1. — : Schematic drawing of the inflated spherical bubble. The spherical bubble is 
powered by jittering jets, i.e., they change their direction at a high rate, launched from an 
accretion disk around the newly formed neutron star. This spherical bubble explodes the 
star according to our model. R-process elements are fused inside the bubble. The typical 
radius during the jets' active phase is 3000-10000 km. 
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thermal energy of the gas and the radiation energy. We neglect losses by neutrinos (see 
Paper 1) and energy production and sink from nuclear reactions. 

During the active time period of the jets the solution to equations (j2J) - (|IJ) is a self- 
similar solution, which we take in the form 

R s (t) = R t a . (5) 

Using equations (jSJ) - (@J we get the following parameters 



a 



3 D/3+5 (/3 + 3)(/3 + 5) 



3 



+ 5' " 12tt^(2/3 + 7)(/3 + 8) J ' lJ 

A short time of ~ 0.1 s after the jets launching process ceases, energy injection to the bubble 
ends. The time delay comes from the jets' crossing time from the NS to the bubble. At 
that moment the self-similar solution no longer holds, and we need to turn to a numerical 
solution. 

The jets' are assumed to be injected at r ~ 15 km fo r a time period of t s = 1 — 2 s. 



Here we take the velocity of the jets to be 150, 000 km s _1 (jCameronll200ll ). This velocity is 
larger than the velocity used in Paper 1 as it better fits the escape velocity from the surface 
of the newly formed NS. The velocity is chosen to be the escape velocity from the neutron 
star, as we are interested in the jets emerging from the neutron star vicinity. Some studies 



inject the jets at distances of > 1000 km at l o wer y elocities, e.g ICouch et al.l (120111 ). while 



others use relativistic jets, e.g., lLazzati et al.l ( 120111 ). who inject the jets at 10 km within 



several seconds. The total mass carried by the two jets is either 0.006M Q or O.OO36M , 
corresponding to a total injected energy of E = 1.7 x 10 51 erg or i? = 1.0 x 10 51 erg 
respectively. For the parameters of E = 1.7 x 10 51 erg s _1 , and active phase time of t s = 1 s, 
for example, the solution during the jets' active phase is 

R 3 (t) = 6.6 x 10 8 t 1 - 3 cm, < t < t s = 1 s. (7) 

For later times we numerically integrate equations (j2J) - (jlj) with Ej = 0, and using the 
the results of the self similar solution at t — t s as initial conditions. For a check, we also 
numerically integrate the equations for the full time of the solution. The numerical solution 
coincides with the self-similar solution after a very short time. The full numerical solution 
for three cases are shown in Fig. [2j The plot shows the radius R s , temperature T, density 
p, and entropy s of the bubble as a function of time. The parameters for the three cases are 
given in the figure caption. The jets' power was taken to be a constant for < t < t s , and 
E = for later times. This gives the little bump at t = t s in the graph. 

The two cases differ by jets' power and having active time of t s = 1 s, are similar in 
their general behavior. As well, changing the active phase duration from t s = 1 s to t s = 2 s 
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does not make large differences in the dynamical properties (the extra density line in the 
left panel of fig. [2j). All cases lead to explosion. Later we will show that these cases are 
different in nucleosynthesis outcomes. 




Fig. 2. — : Left: Radius R s , Temperature T, density p, and entropy s of the spherical bubble 
as a function of time for the model with jets' power of E = 1.7 x 10 51 erg s _1 and a jets' 
active phase lasting t s = 1 s. p\ is the density for a case with Ej = 1 X 10 51 erg s -1 for the 
same duration of T s = 1 s. Right: the same but for a model with E = 0.85 x 10 51 erg s _1 
and an active phase duration of t s = 2 s 

For the typical parameters expected in the model the main results of this section are 
as follows. (1) The spherical bubble reaches a typical radius of ~ 10 4 km at the end of the 
jets' active phase. (2) The temperature relevant for nucleosynthesis (T~3 x 10 9 K) occurs 
at about one seconds from the beginning of the jet injection. (3) The density in the bubble 
of ~ 10 4 g cm" 3 at that time implies that nuclear reactions will be of a high enough rate 
to be significant. For that, in the next section we study the nucleosynthesis inside the bubble. 
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NUCLEOSYNTHESIS INSIDE THE BUBBLE 



In the previous section we found the temperature inside the bubble to start at ~ 10 10 K, 
and to decrease due to adiabatic cooling to 2.5 x 10 9 K in ~ 1 s. The relevant nuclear re- 
actions inside the hot bubble start when the temperature is about T ~ 9 x 10 9 K and 
stop at T ~ 2.5 x 10 9 K. During this time fresh material is injected into the bubble from 
the jets. As the jets are launched fr om very close to t he neutron star, they are composed 
of highly enriched neutron material ( jKohri et al.l 120051 ) . During the expansion of the jets 
they adiabatically cool, and nucleons m i ght fuse to give a par ticles and heavier nuclei (e.g., 
Cameron! 1200 ll ; iMaeda fc Nomotol 120031 ; iFujimoto et al.ll2008l ). However, the jets are even- 
tually shocked with a post shock temperature of > 10 10 K. At that temperature all nuclei 
rapidly disintegrate, and a gas composed of free nucleons is formed. Our nucleosynthesis 
calculations start from the free-nucleons post shock jets' gas. At t ~ 0.2 s when the tem- 
perature inside the spherical bubble has dropped to T ~ 9 x 10 9 K, the free nucleons fusion 
rate overcomes the disintegration rate and a particles start to be accumulated. During the 
time up to ~ 1 s the temperature drops further and a particle fuse to form heavier nuclei 
until a freeze-out is reached. 



The nuclear reaction network is similar to that given in IWooslev fc H offman ( 19921 ). 
The reaction rates are taken from the JINA Reaclib Database ( jCyburt et al.ll2010l ). and 
include reactions with 1,2, and 3 body interactions and beta decays. The reaction network 
is integrated assuming a uniform composition in the bubble and a continues injection of 
protons and neutrons from the jets until the jets terminated. For the ele ctron fraction Y K w e 



<Cohri et al.l d2005h 
Kohri et al.l d2005h 



use values for the accretion disk around a neutron star as calculated by 
For parameters relevant to our model we find from the calculations of 
that the neutron to proton ratio is in the range n/p ~ 5 — 10, namely, Y e ~ 0.09 — 0.17. Here 
we integrate the reaction network for 3 different electron fractions Y e = 0.09, 0.17, 0.25. This 
correspond to a neutron to proton ratio of n/p = 10, 5, and 3. We study nucleosynthesis for 
jets' active phase of 1 s and 2 s. The network is solved independently of the hydrodynamics 
solution of section [5] (IHix fc Thielemann! Il999l ) . 



The nucleosynthesis results are summarized in Fig. |3j The plots show the mass fraction 
of neutrons ,a particles, and of total seed elements during the evolution of the bubble. By 
seed elements we refer to nuclei on which further neutron capt ure will occur to synthesis the 
r-process nuclei ( IWooslev Hoffman! Il992l ; IWitti et al.l I1994J ) . The inclusion of all nuclear 
processes beyond the seed nuclei is beyond the scope of the present paper. They will be 
studied in a forthcoming paper where a full multi-dimensional gasdynamical code will be 
used to study the interaction of the jets with the core material. 



We assume that the post-shock freshly injected jets' material is fully mixed inside the 
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bubble. This is based on the expected formation of vortices inside the bubble by the jittering 
jets. The post-shock velocity for 7 = 4/3 is ~ 20,000 km s -1 , and for a bubble's radius of 
R s ~ 7000 km the mixing time is ~ 0.3 s. This shows that the assumption is reasonable, 
but that mixing is not complete. The assumption of full mixing will be relaxed in the future 
with full multi-dimensional numerical simulations. Since the post-shock entropy increases 
as the jets are shocked at larger distances (because the density is lower), the entropy inside 
the bubble increase as long as we inject fresh jets' material. At the relevant time of nucle- 
osynthesis the entropy per nuc leon reaches values of > 100 Ke/nucleon as is required for 
r-process elements production (jHoffman et al.l 119971 ). Our flow structure differs from calcu- 
latio ns where there is no mixing and the entropy is almost constant during nucleosy nthesis 



e.g. 



Witti et al.lll994! ; IWooslev et al.lll994l lArcones et all 120071 : Kuroda et all 120081 ). 



As well, the continuous injection of nucleons as nucleosynthesis takes place reduces the 
final mass of seed nuclei. This is seen by comparing the mass fraction of seed nuclei for the 
two simulated cases, of 1 s and 2 s jets active phase duration (Table [1] and Fig. [3]). The 
table shows also the sensitivity of the nucleosynthesis production to the value Y e assumed 
at the base of the jet. 







Model I 






Model II 






Model III 




Total mass 


0.60 


0.60 


0.60 


0.36 


0.36 


0.36 


0.60 


0.60 


0.60 


Power 


1.7 


1.7 


1.7 


1.0 


1.0 


1.0 


0.85 


0.85 


0.85 


Ye 


0.09 


0.17 


0.25 


0.09 


0.17 


0.25 


0.09 


0.17 


0.25 


n 


79% 


60% 


42% 


79% 


60% 


42% 


81% 


63% 


46% 


a 


8% 


14% 


21% 


8% 


14% 


20% 


14% 


24% 


33% 


Seed 


13% 


25% 


36% 


13% 


25% 


37% 


4% 


10% 


17% 



Table 1:: Mass fraction of neutrons, a particles and seed elements for different electron 
fractions Y e . Models I and II are for jets with a duration of t s = 1 s. Model III is for jets 
with a duration of t s = 2 s. Total mass is in units of 1O" 2 M . Power is in units of 10 51 erg. 



The most likely case is for jets' actiye phase of 1 — 2 seconds , although the jets' power 



might vary with time (e.g. iHoflich et al.ll200ll ; ICouch et al.ll2009l ). The electron fraction is 



likely to be in the lower part of the range Y P = 0.09 — 0.17, corresponding to a neutron to 
proton ratio of n/p = 5 — 10 (IKohri et al.l 120051 ) . Also, the total energy might be slightly 
less that 1.7 x 10 51 erg used here. From those values we find the seed nuclei mass fraction 
to be 0.02 — 0.2, and the corresponding total seed nuclei mass to be 10 -4 — 10~ 3 M Q , with 
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Fig. 3. — : Evolution of the mass fraction of neutrons, a particles, and seed nuclei for the 
r-process. Left: Model I with t s — 1 s and power of E = 1.7 x 10 51 erg. Right: Model III 
with t s = 2 s and power of E = 0.85 x 10 51 erg. In each panel the initial ratio of electron to 
nucleon Y e is given. The total mass is O.OO6M for both cases. 



more likely values in the range 10~ 4 — 3 x 1O~ 4 M . After neutrons are absorbed (beyond 
the scope of this paper) and form the r-process elements, the total mass of the r-process 
elements is 2 — 3 times that of the seed nuclei. This gives that the expected mass of the 
r-process elements in our jittering-jets explosion model is ~ several x 1O _4 M . 



L 1 



From the solar abundance iMathews &: Cowan! (119901 ) deduced that the average mass 



of r-process material ejected in a CCSN is ~ 10 M . Our simple and idealized model 
overproduces r-process elements, but not by much. In future 3D numerical simulations three 
assumptions that have been used here will be relaxed. These might reduce the production of 
r-process elements by a factor of ~ 3 — 5. (1) Simulating precessing jets will cause deviation 
from sphericity. We expect that in some regions the production of r-process elements will 
be less efficient. (2) Adding the mass at the termination shocks of the jets will result in an 
inhomogeneous bubble. Basically, the flow will differ by having regions where the matter 
cools adiabatically with no addition of entropy (low entropy regions), and regions of high 
entropy where gas is added. This might lead to less efficient production of r-process elements 
in some regions. (3) We will change the the neutron enrichment (or Y e ) of the jets with time. 
It is quite possible that at early times the jets are less neutron enriched than at later time 
when the accretion disk is depleted and more mass comes from closer to the NS. 
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SUMMARY 



Our main goal was to examine the nucleosynthesis inside the bubbles formed by the 
jets in the the jittering- jets model for core collapse SN explosion (Paper 1). The two jets 
are launched from an unrelaxed accretion disk around the newly formed NS. Because of 
stochastic accretion of mass and angular momentum the disk's axis is rapidly changing and 
the disk might even be intermittent. The jets penetrate to a distance of few x 1000 km 
through the infalling stellar core. However, because of their changing direction they cannot 
penetrate beyond that distance. The jets are shocked and form hot low-density bubbles. 
These high pressure bubbles explode the star. This process where the non-penetrating jets 
prevent further accretion to the center (a negative feedback) is termed the penetrating jet 
feedback mechanism. 

To facilitate a solution in the scope of the present paper we assumed that the bubbles 
formed by the two jittering jets merge to form one large spherical bubble, as shown schemat- 
ically in Fig. [TJ The spherical solution under these assumptions is composed of two phases. 
In the first one, the jets' active phase, energy and mass are injected into the bubble at con- 
stant rates. The second phase starts when the jets' cease, and the bubble starts to expand 
adiabatically. The gasdynamical equations in spherical symmetry were solved analytically 
using a self-similar solution for the jets' active phase, and numerically thereafter (section [2]) . 

Beyond the assumptions of jittering jets that have the power to explode the star and 
the formation of a spherical bubble, all quantitative parameters have been used before by 
some studies. We did not adjust or played with any quantitative parameter in the solutions 
presented here. (1) The jets' total energy of 1 — 1. 7xl0 51 erg is taken from the energy required 
to explode CCSNe. (2) The jets' velocity of Vj = 0.5c comes from the escape vel ocity near 
the N S surface. This value for jets' velocity from NS has been used before, e.g., ICameron 
( 120011 ) . (3) The energy and velocity determine the total mass carried by the jets. (4) The 
~ 1 — 2 seco nds duration of the j ets' active phase is similar to durations assumed in other 
studies, e.g. Hoflich et al. ( 2001 h_ Couch et al.l (2009 ). (5) Th e ambient density profile (eq. 
[TJ is taken from lwilson et al" ( 19861 ) and Mikami et al. ( 2008 ). (6) The neutron frac t ion (o r 
electron fraction Y e ) of the jets' material is taken from the calculation of iKohri et al.l ( 120051 ). 



The assumption of a spherical bubble that has the energy to explode the star and the 
quantitative parameters listed above determine the properties and evolution of the bubble. 
From these we calculated the nucleosynthesis inside the bubble. In the limited scope of the 
present paper we numerically integrated a reaction network that follows the fusion up to the 
seed nuclei with Z < 50. The results of the nucleosynthesis calculations for the three studied 
cases are presented in Table [T] and Fig. El We note that Ye = 0:25 is above the expected 
value during the main phase of the jets (IKohri et al.ll2005l ). but might be applicable at early 
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time before the NS is fully relaxed by neutrino cooling. 

During the integration of the network an a freeze-out is reached, i.e., when the number 
of a particles does not change anymore. We take the heavy nuclei from the a freeze-out 
to be the seed elements for r-process elements. From the mass of seed elements for the 
typical parameters expected in this study, Y e = 0.1, energy of 10 51 erg, and active jets' 
duration of t s ~ 1 — 2 s, we estimate the total mass of the fused r-process elements to be 
several x 1O~ 4 M . This is a few times larger th an ~ 10 4 M M , the average mass of r-process 
elements per CCSN deduced from observations ( Mathews &: Cowarilll99ol ). We note that in 



many observed cas es the production of r-process elements is much below the average (e.g. 



Sneden et al.ll2010l and references therein). It is possible, therefore, that the conditions used 
here are met only in a fraction of CCSNe. For example, in some cases the value of Y e at the 
base of the jets is larger than used here, namely Y e > 0.3, where the number of neutrons is 
not sufficient to produce r-process elements. 

One strong assumption of the present work is that the bubble is homogeneous in tem- 
perature and composition. This assumption will be relaxed in a future study when the gas- 
dynamical equations will be solved with a multi-dimensional numerical code. The accurate 
treatment of the inflation process of the bubble will justify the inclusion of a more extended 
nuclear reaction network. Nevertheless, our results here strongly suggest that within the 
context of the jittering-jets model for core collapse SN explosions, the nucleosynthesis of 
the r-process elements is a a likely possibility. This outcome strengthen the possibility that 
CCSNe are driven by jets. 

We thank an anonymous referee for helpful comments. This research was supported by 
the Asher Fund for Space Research at the Technion and the Israel Science foundation. 



REFERENCES 

Arcones, A., Janka, H.-T., & Scheck, L. 2007, A&A, 467, 1227 
Arnould, M., Goriely, S., & Takahashi, K. 2007, Phys. Rep., 450, 97 
Bethe, H. A. 1990, Reviews of Modern Physics, 62, 801 

Bisnovatyi-Kogan, G. S., Popov, I. P., & Samokhin, A. A. 1976, Ap&SS, 41, 287 

Blondin, J. M., & Mezzacappa, A. 2007, Nature, 445, 58 

Brandt, T. D., Burrows, A., Ott, C. D., & Livne, E. 2011, ApJ, 728, 8 



- 12 - 



Cameron, A. G. W. 2001, ApJ, 562, 456 

Chugai, N. N., Fabrika, S. N., Sholukhova, O. N., et al. 2005, Astronomy Letters, 31, 792 
Couch, S. M., Wheeler, J. C, & Milosavljevic, M. 2009, ApJ, 696, 953 
Couch, S. M., Pooley, D., Wheeler, J. C, & Milosavljevic, M. 2011, ApJ, 727, 104 
Cyburt, R. H., et al. 2010, ApJS, 189, 240 

Elmhamdi, A., Danziger, I. J., Chugai, N., et al. 2003, MNRAS, 338, 939 

Fargion, D. 2003, Chinese Journal of Astronomy and Astrophysics Supplement, 3, 472 

Fujimoto, S.-L, Nishimura, N., & Hashimoto, M.-a. 2008, ApJ, 680, 1350 

Hanke, F., Marek, A., Mueller, B., & Janka, H.-T. 2011, arXiv: 1108.43551 

Hix, W. R., & Thielemann, F.-K. 1999, Journal of Computational and Applied Mathematics, 
109, 321 

Hoffman, R. D., Woosley, S. E., & Qian, Y.-Z. 1997, ApJ, 482, 951 

Hoflich, P., Khokhlov, A., & Wang, L. 2001, 20th Texas Symposium on relativistic astro- 
physics, 586, 459 

Khokhlov, A. M., Hdfiich, P. A., Oran, E. S., et al. 1999, ApJ, 524, L107 
Kohri, K., Narayan, R., & Piran, T. 2005, ApJ, 629, 341 
Kuroda, T., Wanajo, S., & Nomoto, K. 2008, ApJ, 672, 1068 

Lazzati, D., Morsony, B. J., Blackwell, C. H., & Begelman, M. C. 2011, arX iv:1111.0970l 
LeBlanc, J. M., & Wilson, J. R. 1970, ApJ, 161, 541 

Leonard, D. C, Filippenko, A. V., Ardila, D. R., & Brotherton, M. S. 2001, ApJ, 553, 861 
Leonard, D. C, Filippenko, A. V., Ganeshalingam, M., et al. 2006, Nature, 440, 505 
MacFadyen, A. I., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 410 
Maeda, K., & Nomoto, K. 2003, ApJ, 598, 1163 
Mathews, G. J., & Cowan, J. J. 1990, Nature, 345, 491 

Meier, D. L., Epstein, R. I., Arnett, W. D., & Schramm, D. N. 1976, ApJ, 204, 869 



-13 - 

Mikami, H., Sato, Y., Matsumoto, T., & Hanawa, T. 2008, ApJ, 683, 357 

Nishimura, S., Kotake, K., Hashimoto, M.-a., Yamada, S., Nishimura, N., Fujimoto, S., & 
Sato, K. 2006, ApJ, 642, 410 

Nordhaus, J., Burrows, A., Almgren, A., & Bell, J. 2010, ApJ, 720, 694 

Papish, O., & Soker, N. 2011, MNRAS, 1321 

Sneden, C, Roederer, I., Cowan, J., & Lawler, J. E. 2010, Nuclei in the Cosmos., 
Smith, N., Cenko, S. B., Butler, N., et al. 2011. larXiv: 1108.28681 
Soker, N. 2010, MNRAS, 401, 2793 
Volk, K., & Kwok, S. 1985, A&A, 153, 79 

Wang, L., Howell, D. A., Hdflich, P., & Wheeler, J. C. 2001, ApJ, 550, 1030 

Wilson, J. R., Mayle, R., Woosley, S. E., & Weaver, T. 1986, Annals of the New York 
Academy of Sciences , 470, 267 

Witti, J., Janka, H.-T., & Takahashi, K. 1994, A&A, 286, 841 

Woosley, S. E., & Hoffman, R. D. 1992, ApJ, 395, 202 

Woosley, S. E., Wilson, J. R., Mathews, G. J., Hoffman, R. D., & Meyer, B. S. 1994, ApJ, 
433, 229 

Woosley, S., & Janka, T. 2005, Nature Physics, 1, 147 



This preprint was prepared with the A AS IATgX macros v5.2. 



